Data preparation and linear trend and monthly means model:
rawd1 <- read.csv("https://montana.box.com/shared/static/s2iudylpfa8gadsnezm7utov2tu9br6l.csv", header= T)
# My first simple idea - generate floor truncated years and then subtract
rawd1$Year <- (floor(rawd1$DATE/100))
rawd1$Month <- round((rawd1$DATE/100 - rawd1$Year) * 100, 1)
rawd1$Yearfrac <- rawd1$Year + (rawd1$Month - 1) / 12
rawd1$Monthf <- factor(rawd1$Month)
model1 <- lm(MMXT ~ Yearfrac + Monthf,data = rawd1)
source("https://montana.box.com/shared/static/u75xpj5lqbp12xuj18vlsn2krrldr9s7.r")
#Function to clean up lm model summaries
print(summary(model1), concise = T)
##
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -69.290162 7.266158 -9.536 <1e-04 ***
## Yearfrac 0.051674 0.003706 13.942 <1e-04 ***
## Monthf2 3.701446 0.604984 6.118 <1e-04 ***
## Monthf3 11.189381 0.604985 18.495 <1e-04 ***
## Monthf4 21.977317 0.604986 36.327 <1e-04 ***
## Monthf5 31.224003 0.606297 51.500 <1e-04 ***
## Monthf6 39.708393 0.606299 65.493 <1e-04 ***
## Monthf7 49.564963 0.608994 81.388 <1e-04 ***
## Monthf8 48.465835 0.607645 79.760 <1e-04 ***
## Monthf9 37.629365 0.608976 61.791 <1e-04 ***
## Monthf10 25.822126 0.607627 42.497 <1e-04 ***
## Monthf11 10.315109 0.607657 16.975 <1e-04 ***
## Monthf12 1.732846 0.608989 2.845 0.0045 **
##
## Residual standard error: 4.597 on 1361 degrees of freedom
## Multiple R-squared: 0.9349, Adjusted R-squared: 0.9343
## F-statistic: 1628 on 12 and 1361 DF, p-value: < 1e-04
# Or using pander:
require(pander)
pander(summary(model1))
| Yearfrac |
0.05167 |
0.003706 |
13.94 |
2.148e-41 |
| Monthf2 |
3.701 |
0.605 |
6.118 |
1.235e-09 |
| Monthf3 |
11.19 |
0.605 |
18.5 |
2.617e-68 |
| Monthf4 |
21.98 |
0.605 |
36.33 |
1.451e-202 |
| Monthf5 |
31.22 |
0.6063 |
51.5 |
6.917e-322 |
| Monthf6 |
39.71 |
0.6063 |
65.49 |
0 |
| Monthf7 |
49.56 |
0.609 |
81.39 |
0 |
| Monthf8 |
48.47 |
0.6076 |
79.76 |
0 |
| Monthf9 |
37.63 |
0.609 |
61.79 |
0 |
| Monthf10 |
25.82 |
0.6076 |
42.5 |
7.23e-252 |
| Monthf11 |
10.32 |
0.6077 |
16.98 |
9.031e-59 |
| Monthf12 |
1.733 |
0.609 |
2.845 |
0.004501 |
| (Intercept) |
-69.29 |
7.266 |
-9.536 |
6.564e-21 |
Fitting linear model: MMXT ~ Yearfrac + Monthf
| 1374 |
4.597 |
0.9349 |
0.9343 |
#There are a few missing months in the data set:
plot(rawd1$Yearfrac,type = "l")
plot(diff(rawd1$Yearfrac) ~ rawd1$Yearfrac[-1],type = "p")
plot(MMXT ~ Yearfrac,data = rawd1,type = "l")
#Complete Yearfrac and Month:
fd <- 1:(115*12 + 11)
ts_complete <- ts(fd,start = 1900,freq = 12)
completetime <- data.frame(Yearfrac = round(as.vector(time(ts_complete)), 4),Month = as.vector(cycle(ts_complete)))
summary(completetime)
## Yearfrac Month
## Min. :1900 Min. : 1.000
## 1st Qu.:1929 1st Qu.: 3.500
## Median :1958 Median : 6.000
## Mean :1958 Mean : 6.496
## 3rd Qu.:1987 3rd Qu.: 9.000
## Max. :2016 Max. :12.000
tail(completetime)
## Yearfrac Month
## 1386 2015.417 6
## 1387 2015.500 7
## 1388 2015.583 8
## 1389 2015.667 9
## 1390 2015.750 10
## 1391 2015.833 11
summary(rawd1)
## STATION STATION_NAME
## COOP:241044:1374 BOZEMAN MONTANA STATE UNIVERSITY MT US:1374
##
##
##
##
##
##
## DATE MMXT Year Month
## Min. :190001 Min. :10.30 Min. :1900 Min. : 1.000
## 1st Qu.:192834 1st Qu.:38.90 1st Qu.:1928 1st Qu.: 3.000
## Median :195806 Median :55.45 Median :1958 Median : 6.000
## Mean :195764 Mean :55.27 Mean :1958 Mean : 6.473
## 3rd Qu.:198679 3rd Qu.:70.60 3rd Qu.:1987 3rd Qu.: 9.000
## Max. :201511 Max. :92.30 Max. :2015 Max. :12.000
##
## Yearfrac Monthf
## Min. :1900 2 :116
## 1st Qu.:1929 3 :116
## Median :1958 4 :116
## Mean :1958 1 :115
## 3rd Qu.:1987 5 :115
## Max. :2016 6 :115
## (Other):681
Fit GAM trend
require(mgcv)
gam1<- gam(MMXT ~ s(Yearfrac,bs = "ts",k = 53) + s(Month, bs ="cc", k = 12), data = rawd1)
gam.check(gam1)
##
## Method: GCV Optimizer: magic
## Smoothing parameter selection converged after 7 iterations.
## The RMS GCV score gradient at convergence was 1.584338e-05 .
## The Hessian was positive definite.
## Model rank = 63 / 63
##
## Basis dimension (k) checking results. Low p-value (k-index<1) may
## indicate that k is too low, especially if edf is close to k'.
##
## k' edf k-index p-value
## s(Yearfrac) 52.00 8.66 0.89 <2e-16 ***
## s(Month) 10.00 9.13 0.99 0.37
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summary(gam1)
##
## Family: gaussian
## Link function: identity
##
## Formula:
## MMXT ~ s(Yearfrac, bs = "ts", k = 53) + s(Month, bs = "cc", k = 12)
##
## Parametric coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 55.2715 0.1236 447 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Approximate significance of smooth terms:
## edf Ref.df F p-value
## s(Yearfrac) 8.660 52 4.131 <2e-16 ***
## s(Month) 9.132 10 1940.301 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## R-sq.(adj) = 0.935 Deviance explained = 93.6%
## GCV = 21.298 Scale est. = 21.006 n = 1374
plot(gam1, scale = -1, pages = 1, shade = T, shade.col = "skyblue")
plot(gam1, scale = 0, pages = 1, shade = T, shade.col = "skyblue")
#Components for sonification at original times of observation:
predres <- predict(gam1, type = "terms")
#Long term trend:
esttrend <- 55.2715 + predres[,1]
plot(esttrend ~ rawd1$Yearfrac, type = "l")
#Seasonal component:
estseasonal <- predres[,2]
plot(estseasonal ~ rawd1$Yearfrac, type = "l") #gap in 50s from missing data...
#Residuals from trend+seasonal model:
resids1 <- residuals(gam1)
plot(resids1 ~ rawd1$Yearfrac, type = "l") #The random unexplained variation in the mean monthly maximum temperatures
par(mfrow = c(3,1))
plot(esttrend ~ rawd1$Yearfrac, type = "l")
plot(estseasonal ~ rawd1$Yearfrac, type = "l")
plot(resids1 ~ rawd1$Yearfrac, type = "l")
#Repeat but make on complete times without gaps:
residsbytime <- data.frame(Yearfrac_obs = round(rawd1$Yearfrac, 4), resids1)
merge1 <- merge(x = completetime, y = residsbytime, by.x = "Yearfrac",by.y = "Yearfrac_obs")
residseries <- matrix(nrow = length(completetime[,1]), ncol = 1,0)
residseries[match(x = residsbytime$Yearfrac,table = completetime$Yearfrac)] <- residsbytime$resids1
par(mfrow = c(1,1))
plot(residseries ~ completetime$Yearfrac, col = as.numeric(residseries == 0) + 1)
#Shows that missing responses were given 0 residuals
#Now to build a complete data set with components:
predres2 <- predict(gam1,newdata = completetime, type = "terms")
#head(predres2)
#head(completecomponents)
esttrend <- 55.2715 + predres2[,1]
estseasonal <- predres2[,2]
completecomponents <- data.frame(esttrend, estseasonal, residseries, Yearfrac = completetime$Yearfrac)
write.csv(completecomponents, "BozemanTemp.csv")
par(mfrow = c(3,1))
Use audiolyzR
library(audiolyzR)
## Loading required package: hexbin
## Loading required package: RJSONIO
## Loading required package: plotrix
audioScatter(esttrend ~ Yearfrac, data = completecomponents)
## R would like permission to write audioplot data files to:
## /Users/huafengzhang during this session.
## Otherwise you will have to manually drag that folder to the synthesizer.
## Is that OK?
## Type y or n
## To load plots to synthesizer:
## Drag the folder entitled "json_matrix"
## to the "Drop folder here to Audiolyze!" box
## at the bottom of the synthesizer (you may have to scroll down).
## Note: Audio plot files written to
## /var/folders/t1/mtpgns5d15n54c8vkk8xm14m0000gn/T//RtmpGyoQfh/json_matrix
audioScatter(estseasonal ~ Yearfrac, data = completecomponents)
## To load plots to synthesizer:
## Drag the folder entitled "json_matrix"
## to the "Drop folder here to Audiolyze!" box
## at the bottom of the synthesizer (you may have to scroll down).
## Note: Audio plot files written to
## /var/folders/t1/mtpgns5d15n54c8vkk8xm14m0000gn/T//RtmpGyoQfh/json_matrix
audioScatter(residseries ~ Yearfrac, data = completecomponents)
## To load plots to synthesizer:
## Drag the folder entitled "json_matrix"
## to the "Drop folder here to Audiolyze!" box
## at the bottom of the synthesizer (you may have to scroll down).
## Note: Audio plot files written to
## /var/folders/t1/mtpgns5d15n54c8vkk8xm14m0000gn/T//RtmpGyoQfh/json_matrix
audioScatter(MMXT ~ Year,data = rawd1)
## To load plots to synthesizer:
## Drag the folder entitled "json_matrix"
## to the "Drop folder here to Audiolyze!" box
## at the bottom of the synthesizer (you may have to scroll down).
## Note: Audio plot files written to
## /var/folders/t1/mtpgns5d15n54c8vkk8xm14m0000gn/T//RtmpGyoQfh/json_matrix
audioScatter(MMXT ~ Month,data = rawd1)
## To load plots to synthesizer:
## Drag the folder entitled "json_matrix"
## to the "Drop folder here to Audiolyze!" box
## at the bottom of the synthesizer (you may have to scroll down).
## Note: Audio plot files written to
## /var/folders/t1/mtpgns5d15n54c8vkk8xm14m0000gn/T//RtmpGyoQfh/json_matrix
audioScatter(MMXT ~ Yearfrac,data = rawd1)
## To load plots to synthesizer:
## Drag the folder entitled "json_matrix"
## to the "Drop folder here to Audiolyze!" box
## at the bottom of the synthesizer (you may have to scroll down).
## Note: Audio plot files written to
## /var/folders/t1/mtpgns5d15n54c8vkk8xm14m0000gn/T//RtmpGyoQfh/json_matrix
Use sonify
library(sonify)
## Loading required package: tuneR
sonify(x = completecomponents$Yearfrac, y = completecomponents$esttrend, waveform ="sine", interpolation ="spline",
duration = 5, flim = c(440, 880), ticks = NULL, tick_len = 0.05,
pulse_len = 0, pulse_amp = 0.2, noise_interval = c(0, 0),
noise_amp = 0.5, amp_level = 1, na_freq = 300, stereo = TRUE,
smp_rate = 44100, play = TRUE, player = NULL, player_args = NULL)
sonify(x = completecomponents$Yearfrac, y = completecomponents$estseasonal, waveform ="sine", interpolation ="spline",
duration = 5, flim = c(440, 880), ticks = NULL, tick_len = 0.05,
pulse_len = 0, pulse_amp = 0.2, noise_interval = c(0, 0),
noise_amp = 0.5, amp_level = 1, na_freq = 300, stereo = TRUE,
smp_rate = 44100, play = TRUE, player = NULL, player_args = NULL)
sonify(x = completecomponents$Yearfrac, y = completecomponents$residseries, waveform ="sine", interpolation ="spline",
duration = 5, flim = c(440, 880), ticks = NULL, tick_len = 0.05,
pulse_len = 0, pulse_amp = 0.2, noise_interval = c(0, 0),
noise_amp = 0.5, amp_level = 1, na_freq = 300, stereo = TRUE,
smp_rate = 44100, play = TRUE, player = NULL, player_args = NULL)
Use seewave
library(seewave)
##
## Attaching package: 'seewave'
## The following object is masked from 'package:plotrix':
##
## rescale
Use audio
library(audio)
##
## Attaching package: 'audio'
## The following object is masked from 'package:tuneR':
##
## play
play_sine <- function(freq, seconds) {
sample_freq <- 4 * freq
samples <- sample_freq * seconds
play(sin(1:samples), sample_freq)
}
temp2freq <- function(temp) {
1000 + 100 * temp
}
play_short <- function(freq) {
a <- play_sine(freq, .1)
wait(a)
}
lapply(temp2freq(rawd1$MMXT), play_short)
## [[1]]
## [1] 1
##
## [[2]]
## [1] 1
##
## [[3]]
## [1] 1
##
## [[4]]
## [1] 1
##
## [[5]]
## [1] 1
##
## [[6]]
## [1] 1
##
## [[7]]
## [1] 1
##
## [[8]]
## [1] 1
##
## [[9]]
## [1] 1
##
## [[10]]
## [1] 1
##
## [[11]]
## [1] 1
##
## [[12]]
## [1] 1
##
## [[13]]
## [1] 1
##
## [[14]]
## [1] 1
##
## [[15]]
## [1] 1
##
## [[16]]
## [1] 1
##
## [[17]]
## [1] 1
##
## [[18]]
## [1] 1
##
## [[19]]
## [1] 1
##
## [[20]]
## [1] 1
##
## [[21]]
## [1] 1
##
## [[22]]
## [1] 1
##
## [[23]]
## [1] 1
##
## [[24]]
## [1] 1
##
## [[25]]
## [1] 1
##
## [[26]]
## [1] 1
##
## [[27]]
## [1] 1
##
## [[28]]
## [1] 1
##
## [[29]]
## [1] 1
##
## [[30]]
## [1] 1
##
## [[31]]
## [1] 1
##
## [[32]]
## [1] 1
##
## [[33]]
## [1] 1
##
## [[34]]
## [1] 1
##
## [[35]]
## [1] 1
##
## [[36]]
## [1] 1
##
## [[37]]
## [1] 1
##
## [[38]]
## [1] 1
##
## [[39]]
## [1] 1
##
## [[40]]
## [1] 1
##
## [[41]]
## [1] 1
##
## [[42]]
## [1] 1
##
## [[43]]
## [1] 1
##
## [[44]]
## [1] 1
##
## [[45]]
## [1] 1
##
## [[46]]
## [1] 1
##
## [[47]]
## [1] 1
##
## [[48]]
## [1] 1
##
## [[49]]
## [1] 1
##
## [[50]]
## [1] 1
##
## [[51]]
## [1] 1
##
## [[52]]
## [1] 1
##
## [[53]]
## [1] 1
##
## [[54]]
## [1] 1
##
## [[55]]
## [1] 1
##
## [[56]]
## [1] 1
##
## [[57]]
## [1] 1
##
## [[58]]
## [1] 1
##
## [[59]]
## [1] 1
##
## [[60]]
## [1] 1
##
## [[61]]
## [1] 1
##
## [[62]]
## [1] 1
##
## [[63]]
## [1] 1
##
## [[64]]
## [1] 1
##
## [[65]]
## [1] 1
##
## [[66]]
## [1] 1
##
## [[67]]
## [1] 1
##
## [[68]]
## [1] 1
##
## [[69]]
## [1] 1
##
## [[70]]
## [1] 1
##
## [[71]]
## [1] 1
##
## [[72]]
## [1] 1
##
## [[73]]
## [1] 1
##
## [[74]]
## [1] 1
##
## [[75]]
## [1] 1
##
## [[76]]
## [1] 1
##
## [[77]]
## [1] 1
##
## [[78]]
## [1] 1
##
## [[79]]
## [1] 1
##
## [[80]]
## [1] 1
##
## [[81]]
## [1] 1
##
## [[82]]
## [1] 1
##
## [[83]]
## [1] 1
##
## [[84]]
## [1] 1
##
## [[85]]
## [1] 1
##
## [[86]]
## [1] 1
##
## [[87]]
## [1] 1
##
## [[88]]
## [1] 1
##
## [[89]]
## [1] 1
##
## [[90]]
## [1] 1
##
## [[91]]
## [1] 1
##
## [[92]]
## [1] 1
##
## [[93]]
## [1] 1
##
## [[94]]
## [1] 1
##
## [[95]]
## [1] 1
##
## [[96]]
## [1] 1
##
## [[97]]
## [1] 1
##
## [[98]]
## [1] 1
##
## [[99]]
## [1] 1
##
## [[100]]
## [1] 1
##
## [[101]]
## [1] 1
##
## [[102]]
## [1] 1
##
## [[103]]
## [1] 1
##
## [[104]]
## [1] 1
##
## [[105]]
## [1] 1
##
## [[106]]
## [1] 1
##
## [[107]]
## [1] 1
##
## [[108]]
## [1] 1
##
## [[109]]
## [1] 1
##
## [[110]]
## [1] 1
##
## [[111]]
## [1] 1
##
## [[112]]
## [1] 1
##
## [[113]]
## [1] 1
##
## [[114]]
## [1] 1
##
## [[115]]
## [1] 1
##
## [[116]]
## [1] 1
##
## [[117]]
## [1] 1
##
## [[118]]
## [1] 1
##
## [[119]]
## [1] 1
##
## [[120]]
## [1] 1
##
## [[121]]
## [1] 1
##
## [[122]]
## [1] 1
##
## [[123]]
## [1] 1
##
## [[124]]
## [1] 1
##
## [[125]]
## [1] 1
##
## [[126]]
## [1] 1
##
## [[127]]
## [1] 1
##
## [[128]]
## [1] 1
##
## [[129]]
## [1] 1
##
## [[130]]
## [1] 1
##
## [[131]]
## [1] 1
##
## [[132]]
## [1] 1
##
## [[133]]
## [1] 1
##
## [[134]]
## [1] 1
##
## [[135]]
## [1] 1
##
## [[136]]
## [1] 1
##
## [[137]]
## [1] 1
##
## [[138]]
## [1] 1
##
## [[139]]
## [1] 1
##
## [[140]]
## [1] 1
##
## [[141]]
## [1] 1
##
## [[142]]
## [1] 1
##
## [[143]]
## [1] 1
##
## [[144]]
## [1] 1
##
## [[145]]
## [1] 1
##
## [[146]]
## [1] 1
##
## [[147]]
## [1] 1
##
## [[148]]
## [1] 1
##
## [[149]]
## [1] 1
##
## [[150]]
## [1] 1
##
## [[151]]
## [1] 1
##
## [[152]]
## [1] 1
##
## [[153]]
## [1] 1
##
## [[154]]
## [1] 1
##
## [[155]]
## [1] 1
##
## [[156]]
## [1] 1
##
## [[157]]
## [1] 1
##
## [[158]]
## [1] 1
##
## [[159]]
## [1] 1
##
## [[160]]
## [1] 1
##
## [[161]]
## [1] 1
##
## [[162]]
## [1] 1
##
## [[163]]
## [1] 1
##
## [[164]]
## [1] 1
##
## [[165]]
## [1] 1
##
## [[166]]
## [1] 1
##
## [[167]]
## [1] 1
##
## [[168]]
## [1] 1
##
## [[169]]
## [1] 1
##
## [[170]]
## [1] 1
##
## [[171]]
## [1] 1
##
## [[172]]
## [1] 1
##
## [[173]]
## [1] 1
##
## [[174]]
## [1] 1
##
## [[175]]
## [1] 1
##
## [[176]]
## [1] 1
##
## [[177]]
## [1] 1
##
## [[178]]
## [1] 1
##
## [[179]]
## [1] 1
##
## [[180]]
## [1] 1
##
## [[181]]
## [1] 1
##
## [[182]]
## [1] 1
##
## [[183]]
## [1] 1
##
## [[184]]
## [1] 1
##
## [[185]]
## [1] 1
##
## [[186]]
## [1] 1
##
## [[187]]
## [1] 1
##
## [[188]]
## [1] 1
##
## [[189]]
## [1] 1
##
## [[190]]
## [1] 1
##
## [[191]]
## [1] 1
##
## [[192]]
## [1] 1
##
## [[193]]
## [1] 1
##
## [[194]]
## [1] 1
##
## [[195]]
## [1] 1
##
## [[196]]
## [1] 1
##
## [[197]]
## [1] 1
##
## [[198]]
## [1] 1
##
## [[199]]
## [1] 1
##
## [[200]]
## [1] 1
##
## [[201]]
## [1] 1
##
## [[202]]
## [1] 1
##
## [[203]]
## [1] 1
##
## [[204]]
## [1] 1
##
## [[205]]
## [1] 1
##
## [[206]]
## [1] 1
##
## [[207]]
## [1] 1
##
## [[208]]
## [1] 1
##
## [[209]]
## [1] 1
##
## [[210]]
## [1] 1
##
## [[211]]
## [1] 1
##
## [[212]]
## [1] 1
##
## [[213]]
## [1] 1
##
## [[214]]
## [1] 1
##
## [[215]]
## [1] 1
##
## [[216]]
## [1] 1
##
## [[217]]
## [1] 1
##
## [[218]]
## [1] 1
##
## [[219]]
## [1] 1
##
## [[220]]
## [1] 1
##
## [[221]]
## [1] 1
##
## [[222]]
## [1] 1
##
## [[223]]
## [1] 1
##
## [[224]]
## [1] 1
##
## [[225]]
## [1] 1
##
## [[226]]
## [1] 1
##
## [[227]]
## [1] 1
##
## [[228]]
## [1] 1
##
## [[229]]
## [1] 1
##
## [[230]]
## [1] 1
##
## [[231]]
## [1] 1
##
## [[232]]
## [1] 1
##
## [[233]]
## [1] 1
##
## [[234]]
## [1] 1
##
## [[235]]
## [1] 1
##
## [[236]]
## [1] 1
##
## [[237]]
## [1] 1
##
## [[238]]
## [1] 1
##
## [[239]]
## [1] 1
##
## [[240]]
## [1] 1
##
## [[241]]
## [1] 1
##
## [[242]]
## [1] 1
##
## [[243]]
## [1] 1
##
## [[244]]
## [1] 1
##
## [[245]]
## [1] 1
##
## [[246]]
## [1] 1
##
## [[247]]
## [1] 1
##
## [[248]]
## [1] 1
##
## [[249]]
## [1] 1
##
## [[250]]
## [1] 1
##
## [[251]]
## [1] 1
##
## [[252]]
## [1] 1
##
## [[253]]
## [1] 1
##
## [[254]]
## [1] 1
##
## [[255]]
## [1] 1
##
## [[256]]
## [1] 1
##
## [[257]]
## [1] 1
##
## [[258]]
## [1] 1
##
## [[259]]
## [1] 1
##
## [[260]]
## [1] 1
##
## [[261]]
## [1] 1
##
## [[262]]
## [1] 1
##
## [[263]]
## [1] 1
##
## [[264]]
## [1] 1
##
## [[265]]
## [1] 1
##
## [[266]]
## [1] 1
##
## [[267]]
## [1] 1
##
## [[268]]
## [1] 1
##
## [[269]]
## [1] 1
##
## [[270]]
## [1] 1
##
## [[271]]
## [1] 1
##
## [[272]]
## [1] 1
##
## [[273]]
## [1] 1
##
## [[274]]
## [1] 1
##
## [[275]]
## [1] 1
##
## [[276]]
## [1] 1
##
## [[277]]
## [1] 1
##
## [[278]]
## [1] 1
##
## [[279]]
## [1] 1
##
## [[280]]
## [1] 1
##
## [[281]]
## [1] 1
##
## [[282]]
## [1] 1
##
## [[283]]
## [1] 1
##
## [[284]]
## [1] 1
##
## [[285]]
## [1] 1
##
## [[286]]
## [1] 1
##
## [[287]]
## [1] 1
##
## [[288]]
## [1] 1
##
## [[289]]
## [1] 1
##
## [[290]]
## [1] 1
##
## [[291]]
## [1] 1
##
## [[292]]
## [1] 1
##
## [[293]]
## [1] 1
##
## [[294]]
## [1] 1
##
## [[295]]
## [1] 1
##
## [[296]]
## [1] 1
##
## [[297]]
## [1] 1
##
## [[298]]
## [1] 1
##
## [[299]]
## [1] 1
##
## [[300]]
## [1] 1
##
## [[301]]
## [1] 1
##
## [[302]]
## [1] 1
##
## [[303]]
## [1] 1
##
## [[304]]
## [1] 1
##
## [[305]]
## [1] 1
##
## [[306]]
## [1] 1
##
## [[307]]
## [1] 1
##
## [[308]]
## [1] 1
##
## [[309]]
## [1] 1
##
## [[310]]
## [1] 1
##
## [[311]]
## [1] 1
##
## [[312]]
## [1] 1
##
## [[313]]
## [1] 1
##
## [[314]]
## [1] 1
##
## [[315]]
## [1] 1
##
## [[316]]
## [1] 1
##
## [[317]]
## [1] 1
##
## [[318]]
## [1] 1
##
## [[319]]
## [1] 1
##
## [[320]]
## [1] 1
##
## [[321]]
## [1] 1
##
## [[322]]
## [1] 1
##
## [[323]]
## [1] 1
##
## [[324]]
## [1] 1
##
## [[325]]
## [1] 1
##
## [[326]]
## [1] 1
##
## [[327]]
## [1] 1
##
## [[328]]
## [1] 1
##
## [[329]]
## [1] 1
##
## [[330]]
## [1] 1
##
## [[331]]
## [1] 1
##
## [[332]]
## [1] 1
##
## [[333]]
## [1] 1
##
## [[334]]
## [1] 1
##
## [[335]]
## [1] 1
##
## [[336]]
## [1] 1
##
## [[337]]
## [1] 1
##
## [[338]]
## [1] 1
##
## [[339]]
## [1] 1
##
## [[340]]
## [1] 1
##
## [[341]]
## [1] 1
##
## [[342]]
## [1] 1
##
## [[343]]
## [1] 1
##
## [[344]]
## [1] 1
##
## [[345]]
## [1] 1
##
## [[346]]
## [1] 1
##
## [[347]]
## [1] 1
##
## [[348]]
## [1] 1
##
## [[349]]
## [1] 1
##
## [[350]]
## [1] 1
##
## [[351]]
## [1] 1
##
## [[352]]
## [1] 1
##
## [[353]]
## [1] 1
##
## [[354]]
## [1] 1
##
## [[355]]
## [1] 1
##
## [[356]]
## [1] 1
##
## [[357]]
## [1] 1
##
## [[358]]
## [1] 1
##
## [[359]]
## [1] 1
##
## [[360]]
## [1] 1
##
## [[361]]
## [1] 1
##
## [[362]]
## [1] 1
##
## [[363]]
## [1] 1
##
## [[364]]
## [1] 1
##
## [[365]]
## [1] 1
##
## [[366]]
## [1] 1
##
## [[367]]
## [1] 1
##
## [[368]]
## [1] 1
##
## [[369]]
## [1] 1
##
## [[370]]
## [1] 1
##
## [[371]]
## [1] 1
##
## [[372]]
## [1] 1
##
## [[373]]
## [1] 1
##
## [[374]]
## [1] 1
##
## [[375]]
## [1] 1
##
## [[376]]
## [1] 1
##
## [[377]]
## [1] 1
##
## [[378]]
## [1] 1
##
## [[379]]
## [1] 1
##
## [[380]]
## [1] 1
##
## [[381]]
## [1] 1
##
## [[382]]
## [1] 1
##
## [[383]]
## [1] 1
##
## [[384]]
## [1] 1
##
## [[385]]
## [1] 1
##
## [[386]]
## [1] 1
##
## [[387]]
## [1] 1
##
## [[388]]
## [1] 1
##
## [[389]]
## [1] 1
##
## [[390]]
## [1] 1
##
## [[391]]
## [1] 1
##
## [[392]]
## [1] 1
##
## [[393]]
## [1] 1
##
## [[394]]
## [1] 1
##
## [[395]]
## [1] 1
##
## [[396]]
## [1] 1
##
## [[397]]
## [1] 1
##
## [[398]]
## [1] 1
##
## [[399]]
## [1] 1
##
## [[400]]
## [1] 1
##
## [[401]]
## [1] 1
##
## [[402]]
## [1] 1
##
## [[403]]
## [1] 1
##
## [[404]]
## [1] 1
##
## [[405]]
## [1] 1
##
## [[406]]
## [1] 1
##
## [[407]]
## [1] 1
##
## [[408]]
## [1] 1
##
## [[409]]
## [1] 1
##
## [[410]]
## [1] 1
##
## [[411]]
## [1] 1
##
## [[412]]
## [1] 1
##
## [[413]]
## [1] 1
##
## [[414]]
## [1] 1
##
## [[415]]
## [1] 1
##
## [[416]]
## [1] 1
##
## [[417]]
## [1] 1
##
## [[418]]
## [1] 1
##
## [[419]]
## [1] 1
##
## [[420]]
## [1] 1
##
## [[421]]
## [1] 1
##
## [[422]]
## [1] 1
##
## [[423]]
## [1] 1
##
## [[424]]
## [1] 1
##
## [[425]]
## [1] 1
##
## [[426]]
## [1] 1
##
## [[427]]
## [1] 1
##
## [[428]]
## [1] 1
##
## [[429]]
## [1] 1
##
## [[430]]
## [1] 1
##
## [[431]]
## [1] 1
##
## [[432]]
## [1] 1
##
## [[433]]
## [1] 1
##
## [[434]]
## [1] 1
##
## [[435]]
## [1] 1
##
## [[436]]
## [1] 1
##
## [[437]]
## [1] 1
##
## [[438]]
## [1] 1
##
## [[439]]
## [1] 1
##
## [[440]]
## [1] 1
##
## [[441]]
## [1] 1
##
## [[442]]
## [1] 1
##
## [[443]]
## [1] 1
##
## [[444]]
## [1] 1
##
## [[445]]
## [1] 1
##
## [[446]]
## [1] 1
##
## [[447]]
## [1] 1
##
## [[448]]
## [1] 1
##
## [[449]]
## [1] 1
##
## [[450]]
## [1] 1
##
## [[451]]
## [1] 1
##
## [[452]]
## [1] 1
##
## [[453]]
## [1] 1
##
## [[454]]
## [1] 1
##
## [[455]]
## [1] 1
##
## [[456]]
## [1] 1
##
## [[457]]
## [1] 1
##
## [[458]]
## [1] 1
##
## [[459]]
## [1] 1
##
## [[460]]
## [1] 1
##
## [[461]]
## [1] 1
##
## [[462]]
## [1] 1
##
## [[463]]
## [1] 1
##
## [[464]]
## [1] 1
##
## [[465]]
## [1] 1
##
## [[466]]
## [1] 1
##
## [[467]]
## [1] 1
##
## [[468]]
## [1] 1
##
## [[469]]
## [1] 1
##
## [[470]]
## [1] 1
##
## [[471]]
## [1] 1
##
## [[472]]
## [1] 1
##
## [[473]]
## [1] 1
##
## [[474]]
## [1] 1
##
## [[475]]
## [1] 1
##
## [[476]]
## [1] 1
##
## [[477]]
## [1] 1
##
## [[478]]
## [1] 1
##
## [[479]]
## [1] 1
##
## [[480]]
## [1] 1
##
## [[481]]
## [1] 1
##
## [[482]]
## [1] 1
##
## [[483]]
## [1] 1
##
## [[484]]
## [1] 1
##
## [[485]]
## [1] 1
##
## [[486]]
## [1] 1
##
## [[487]]
## [1] 1
##
## [[488]]
## [1] 1
##
## [[489]]
## [1] 1
##
## [[490]]
## [1] 1
##
## [[491]]
## [1] 1
##
## [[492]]
## [1] 1
##
## [[493]]
## [1] 1
##
## [[494]]
## [1] 1
##
## [[495]]
## [1] 1
##
## [[496]]
## [1] 1
##
## [[497]]
## [1] 1
##
## [[498]]
## [1] 1
##
## [[499]]
## [1] 1
##
## [[500]]
## [1] 1
##
## [[501]]
## [1] 1
##
## [[502]]
## [1] 1
##
## [[503]]
## [1] 1
##
## [[504]]
## [1] 1
##
## [[505]]
## [1] 1
##
## [[506]]
## [1] 1
##
## [[507]]
## [1] 1
##
## [[508]]
## [1] 1
##
## [[509]]
## [1] 1
##
## [[510]]
## [1] 1
##
## [[511]]
## [1] 1
##
## [[512]]
## [1] 1
##
## [[513]]
## [1] 1
##
## [[514]]
## [1] 1
##
## [[515]]
## [1] 1
##
## [[516]]
## [1] 1
##
## [[517]]
## [1] 1
##
## [[518]]
## [1] 1
##
## [[519]]
## [1] 1
##
## [[520]]
## [1] 1
##
## [[521]]
## [1] 1
##
## [[522]]
## [1] 1
##
## [[523]]
## [1] 1
##
## [[524]]
## [1] 1
##
## [[525]]
## [1] 1
##
## [[526]]
## [1] 1
##
## [[527]]
## [1] 1
##
## [[528]]
## [1] 1
##
## [[529]]
## [1] 1
##
## [[530]]
## [1] 1
##
## [[531]]
## [1] 1
##
## [[532]]
## [1] 1
##
## [[533]]
## [1] 1
##
## [[534]]
## [1] 1
##
## [[535]]
## [1] 1
##
## [[536]]
## [1] 1
##
## [[537]]
## [1] 1
##
## [[538]]
## [1] 1
##
## [[539]]
## [1] 1
##
## [[540]]
## [1] 1
##
## [[541]]
## [1] 1
##
## [[542]]
## [1] 1
##
## [[543]]
## [1] 1
##
## [[544]]
## [1] 1
##
## [[545]]
## [1] 1
##
## [[546]]
## [1] 1
##
## [[547]]
## [1] 1
##
## [[548]]
## [1] 1
##
## [[549]]
## [1] 1
##
## [[550]]
## [1] 1
##
## [[551]]
## [1] 1
##
## [[552]]
## [1] 1
##
## [[553]]
## [1] 1
##
## [[554]]
## [1] 1
##
## [[555]]
## [1] 1
##
## [[556]]
## [1] 1
##
## [[557]]
## [1] 1
##
## [[558]]
## [1] 1
##
## [[559]]
## [1] 1
##
## [[560]]
## [1] 1
##
## [[561]]
## [1] 1
##
## [[562]]
## [1] 1
##
## [[563]]
## [1] 1
##
## [[564]]
## [1] 1
##
## [[565]]
## [1] 1
##
## [[566]]
## [1] 1
##
## [[567]]
## [1] 1
##
## [[568]]
## [1] 1
##
## [[569]]
## [1] 1
##
## [[570]]
## [1] 1
##
## [[571]]
## [1] 1
##
## [[572]]
## [1] 1
##
## [[573]]
## [1] 1
##
## [[574]]
## [1] 1
##
## [[575]]
## [1] 1
##
## [[576]]
## [1] 1
##
## [[577]]
## [1] 1
##
## [[578]]
## [1] 1
##
## [[579]]
## [1] 1
##
## [[580]]
## [1] 1
##
## [[581]]
## [1] 1
##
## [[582]]
## [1] 1
##
## [[583]]
## [1] 1
##
## [[584]]
## [1] 1
##
## [[585]]
## [1] 1
##
## [[586]]
## [1] 1
##
## [[587]]
## [1] 1
##
## [[588]]
## [1] 1
##
## [[589]]
## [1] 1
##
## [[590]]
## [1] 1
##
## [[591]]
## [1] 1
##
## [[592]]
## [1] 1
##
## [[593]]
## [1] 1
##
## [[594]]
## [1] 1
##
## [[595]]
## [1] 1
##
## [[596]]
## [1] 1
##
## [[597]]
## [1] 1
##
## [[598]]
## [1] 1
##
## [[599]]
## [1] 1
##
## [[600]]
## [1] 1
##
## [[601]]
## [1] 1
##
## [[602]]
## [1] 1
##
## [[603]]
## [1] 1
##
## [[604]]
## [1] 1
##
## [[605]]
## [1] 1
##
## [[606]]
## [1] 1
##
## [[607]]
## [1] 1
##
## [[608]]
## [1] 1
##
## [[609]]
## [1] 1
##
## [[610]]
## [1] 1
##
## [[611]]
## [1] 1
##
## [[612]]
## [1] 1
##
## [[613]]
## [1] 1
##
## [[614]]
## [1] 1
##
## [[615]]
## [1] 1
##
## [[616]]
## [1] 1
##
## [[617]]
## [1] 1
##
## [[618]]
## [1] 1
##
## [[619]]
## [1] 1
##
## [[620]]
## [1] 1
##
## [[621]]
## [1] 1
##
## [[622]]
## [1] 1
##
## [[623]]
## [1] 1
##
## [[624]]
## [1] 1
##
## [[625]]
## [1] 1
##
## [[626]]
## [1] 1
##
## [[627]]
## [1] 1
##
## [[628]]
## [1] 1
##
## [[629]]
## [1] 1
##
## [[630]]
## [1] 1
##
## [[631]]
## [1] 1
##
## [[632]]
## [1] 1
##
## [[633]]
## [1] 1
##
## [[634]]
## [1] 1
##
## [[635]]
## [1] 1
##
## [[636]]
## [1] 1
##
## [[637]]
## [1] 1
##
## [[638]]
## [1] 1
##
## [[639]]
## [1] 1
##
## [[640]]
## [1] 1
##
## [[641]]
## [1] 1
##
## [[642]]
## [1] 1
##
## [[643]]
## [1] 1
##
## [[644]]
## [1] 1
##
## [[645]]
## [1] 1
##
## [[646]]
## [1] 1
##
## [[647]]
## [1] 1
##
## [[648]]
## [1] 1
##
## [[649]]
## [1] 1
##
## [[650]]
## [1] 1
##
## [[651]]
## [1] 1
##
## [[652]]
## [1] 1
##
## [[653]]
## [1] 1
##
## [[654]]
## [1] 1
##
## [[655]]
## [1] 1
##
## [[656]]
## [1] 1
##
## [[657]]
## [1] 1
##
## [[658]]
## [1] 1
##
## [[659]]
## [1] 1
##
## [[660]]
## [1] 1
##
## [[661]]
## [1] 1
##
## [[662]]
## [1] 1
##
## [[663]]
## [1] 1
##
## [[664]]
## [1] 1
##
## [[665]]
## [1] 1
##
## [[666]]
## [1] 1
##
## [[667]]
## [1] 1
##
## [[668]]
## [1] 1
##
## [[669]]
## [1] 1
##
## [[670]]
## [1] 1
##
## [[671]]
## [1] 1
##
## [[672]]
## [1] 1
##
## [[673]]
## [1] 1
##
## [[674]]
## [1] 1
##
## [[675]]
## [1] 1
##
## [[676]]
## [1] 1
##
## [[677]]
## [1] 1
##
## [[678]]
## [1] 1
##
## [[679]]
## [1] 1
##
## [[680]]
## [1] 1
##
## [[681]]
## [1] 1
##
## [[682]]
## [1] 1
##
## [[683]]
## [1] 1
##
## [[684]]
## [1] 1
##
## [[685]]
## [1] 1
##
## [[686]]
## [1] 1
##
## [[687]]
## [1] 1
##
## [[688]]
## [1] 1
##
## [[689]]
## [1] 1
##
## [[690]]
## [1] 1
##
## [[691]]
## [1] 1
##
## [[692]]
## [1] 1
##
## [[693]]
## [1] 1
##
## [[694]]
## [1] 1
##
## [[695]]
## [1] 1
##
## [[696]]
## [1] 1
##
## [[697]]
## [1] 1
##
## [[698]]
## [1] 1
##
## [[699]]
## [1] 1
##
## [[700]]
## [1] 1
##
## [[701]]
## [1] 1
##
## [[702]]
## [1] 1
##
## [[703]]
## [1] 1
##
## [[704]]
## [1] 1
##
## [[705]]
## [1] 1
##
## [[706]]
## [1] 1
##
## [[707]]
## [1] 1
##
## [[708]]
## [1] 1
##
## [[709]]
## [1] 1
##
## [[710]]
## [1] 1
##
## [[711]]
## [1] 1
##
## [[712]]
## [1] 1
##
## [[713]]
## [1] 1
##
## [[714]]
## [1] 1
##
## [[715]]
## [1] 1
##
## [[716]]
## [1] 1
##
## [[717]]
## [1] 1
##
## [[718]]
## [1] 1
##
## [[719]]
## [1] 1
##
## [[720]]
## [1] 1
##
## [[721]]
## [1] 1
##
## [[722]]
## [1] 1
##
## [[723]]
## [1] 1
##
## [[724]]
## [1] 1
##
## [[725]]
## [1] 1
##
## [[726]]
## [1] 1
##
## [[727]]
## [1] 1
##
## [[728]]
## [1] 1
##
## [[729]]
## [1] 1
##
## [[730]]
## [1] 1
##
## [[731]]
## [1] 1
##
## [[732]]
## [1] 1
##
## [[733]]
## [1] 1
##
## [[734]]
## [1] 1
##
## [[735]]
## [1] 1
##
## [[736]]
## [1] 1
##
## [[737]]
## [1] 1
##
## [[738]]
## [1] 1
##
## [[739]]
## [1] 1
##
## [[740]]
## [1] 1
##
## [[741]]
## [1] 1
##
## [[742]]
## [1] 1
##
## [[743]]
## [1] 1
##
## [[744]]
## [1] 1
##
## [[745]]
## [1] 1
##
## [[746]]
## [1] 1
##
## [[747]]
## [1] 1
##
## [[748]]
## [1] 1
##
## [[749]]
## [1] 1
##
## [[750]]
## [1] 1
##
## [[751]]
## [1] 1
##
## [[752]]
## [1] 1
##
## [[753]]
## [1] 1
##
## [[754]]
## [1] 1
##
## [[755]]
## [1] 1
##
## [[756]]
## [1] 1
##
## [[757]]
## [1] 1
##
## [[758]]
## [1] 1
##
## [[759]]
## [1] 1
##
## [[760]]
## [1] 1
##
## [[761]]
## [1] 1
##
## [[762]]
## [1] 1
##
## [[763]]
## [1] 1
##
## [[764]]
## [1] 1
##
## [[765]]
## [1] 1
##
## [[766]]
## [1] 1
##
## [[767]]
## [1] 1
##
## [[768]]
## [1] 1
##
## [[769]]
## [1] 1
##
## [[770]]
## [1] 1
##
## [[771]]
## [1] 1
##
## [[772]]
## [1] 1
##
## [[773]]
## [1] 1
##
## [[774]]
## [1] 1
##
## [[775]]
## [1] 1
##
## [[776]]
## [1] 1
##
## [[777]]
## [1] 1
##
## [[778]]
## [1] 1
##
## [[779]]
## [1] 1
##
## [[780]]
## [1] 1
##
## [[781]]
## [1] 1
##
## [[782]]
## [1] 1
##
## [[783]]
## [1] 1
##
## [[784]]
## [1] 1
##
## [[785]]
## [1] 1
##
## [[786]]
## [1] 1
##
## [[787]]
## [1] 1
##
## [[788]]
## [1] 1
##
## [[789]]
## [1] 1
##
## [[790]]
## [1] 1
##
## [[791]]
## [1] 1
##
## [[792]]
## [1] 1
##
## [[793]]
## [1] 1
##
## [[794]]
## [1] 1
##
## [[795]]
## [1] 1
##
## [[796]]
## [1] 1
##
## [[797]]
## [1] 1
##
## [[798]]
## [1] 1
##
## [[799]]
## [1] 1
##
## [[800]]
## [1] 1
##
## [[801]]
## [1] 1
##
## [[802]]
## [1] 1
##
## [[803]]
## [1] 1
##
## [[804]]
## [1] 1
##
## [[805]]
## [1] 1
##
## [[806]]
## [1] 1
##
## [[807]]
## [1] 1
##
## [[808]]
## [1] 1
##
## [[809]]
## [1] 1
##
## [[810]]
## [1] 1
##
## [[811]]
## [1] 1
##
## [[812]]
## [1] 1
##
## [[813]]
## [1] 1
##
## [[814]]
## [1] 1
##
## [[815]]
## [1] 1
##
## [[816]]
## [1] 1
##
## [[817]]
## [1] 1
##
## [[818]]
## [1] 1
##
## [[819]]
## [1] 1
##
## [[820]]
## [1] 1
##
## [[821]]
## [1] 1
##
## [[822]]
## [1] 1
##
## [[823]]
## [1] 1
##
## [[824]]
## [1] 1
##
## [[825]]
## [1] 1
##
## [[826]]
## [1] 1
##
## [[827]]
## [1] 1
##
## [[828]]
## [1] 1
##
## [[829]]
## [1] 1
##
## [[830]]
## [1] 1
##
## [[831]]
## [1] 1
##
## [[832]]
## [1] 1
##
## [[833]]
## [1] 1
##
## [[834]]
## [1] 1
##
## [[835]]
## [1] 1
##
## [[836]]
## [1] 1
##
## [[837]]
## [1] 1
##
## [[838]]
## [1] 1
##
## [[839]]
## [1] 1
##
## [[840]]
## [1] 1
##
## [[841]]
## [1] 1
##
## [[842]]
## [1] 1
##
## [[843]]
## [1] 1
##
## [[844]]
## [1] 1
##
## [[845]]
## [1] 1
##
## [[846]]
## [1] 1
##
## [[847]]
## [1] 1
##
## [[848]]
## [1] 1
##
## [[849]]
## [1] 1
##
## [[850]]
## [1] 1
##
## [[851]]
## [1] 1
##
## [[852]]
## [1] 1
##
## [[853]]
## [1] 1
##
## [[854]]
## [1] 1
##
## [[855]]
## [1] 1
##
## [[856]]
## [1] 1
##
## [[857]]
## [1] 1
##
## [[858]]
## [1] 1
##
## [[859]]
## [1] 1
##
## [[860]]
## [1] 1
##
## [[861]]
## [1] 1
##
## [[862]]
## [1] 1
##
## [[863]]
## [1] 1
##
## [[864]]
## [1] 1
##
## [[865]]
## [1] 1
##
## [[866]]
## [1] 1
##
## [[867]]
## [1] 1
##
## [[868]]
## [1] 1
##
## [[869]]
## [1] 1
##
## [[870]]
## [1] 1
##
## [[871]]
## [1] 1
##
## [[872]]
## [1] 1
##
## [[873]]
## [1] 1
##
## [[874]]
## [1] 1
##
## [[875]]
## [1] 1
##
## [[876]]
## [1] 1
##
## [[877]]
## [1] 1
##
## [[878]]
## [1] 1
##
## [[879]]
## [1] 1
##
## [[880]]
## [1] 1
##
## [[881]]
## [1] 1
##
## [[882]]
## [1] 1
##
## [[883]]
## [1] 1
##
## [[884]]
## [1] 1
##
## [[885]]
## [1] 1
##
## [[886]]
## [1] 1
##
## [[887]]
## [1] 1
##
## [[888]]
## [1] 1
##
## [[889]]
## [1] 1
##
## [[890]]
## [1] 1
##
## [[891]]
## [1] 1
##
## [[892]]
## [1] 1
##
## [[893]]
## [1] 1
##
## [[894]]
## [1] 1
##
## [[895]]
## [1] 1
##
## [[896]]
## [1] 1
##
## [[897]]
## [1] 1
##
## [[898]]
## [1] 1
##
## [[899]]
## [1] 1
##
## [[900]]
## [1] 1
##
## [[901]]
## [1] 1
##
## [[902]]
## [1] 1
##
## [[903]]
## [1] 1
##
## [[904]]
## [1] 1
##
## [[905]]
## [1] 1
##
## [[906]]
## [1] 1
##
## [[907]]
## [1] 1
##
## [[908]]
## [1] 1
##
## [[909]]
## [1] 1
##
## [[910]]
## [1] 1
##
## [[911]]
## [1] 1
##
## [[912]]
## [1] 1
##
## [[913]]
## [1] 1
##
## [[914]]
## [1] 1
##
## [[915]]
## [1] 1
##
## [[916]]
## [1] 1
##
## [[917]]
## [1] 1
##
## [[918]]
## [1] 1
##
## [[919]]
## [1] 1
##
## [[920]]
## [1] 1
##
## [[921]]
## [1] 1
##
## [[922]]
## [1] 1
##
## [[923]]
## [1] 1
##
## [[924]]
## [1] 1
##
## [[925]]
## [1] 1
##
## [[926]]
## [1] 1
##
## [[927]]
## [1] 1
##
## [[928]]
## [1] 1
##
## [[929]]
## [1] 1
##
## [[930]]
## [1] 1
##
## [[931]]
## [1] 1
##
## [[932]]
## [1] 1
##
## [[933]]
## [1] 1
##
## [[934]]
## [1] 1
##
## [[935]]
## [1] 1
##
## [[936]]
## [1] 1
##
## [[937]]
## [1] 1
##
## [[938]]
## [1] 1
##
## [[939]]
## [1] 1
##
## [[940]]
## [1] 1
##
## [[941]]
## [1] 1
##
## [[942]]
## [1] 1
##
## [[943]]
## [1] 1
##
## [[944]]
## [1] 1
##
## [[945]]
## [1] 1
##
## [[946]]
## [1] 1
##
## [[947]]
## [1] 1
##
## [[948]]
## [1] 1
##
## [[949]]
## [1] 1
##
## [[950]]
## [1] 1
##
## [[951]]
## [1] 1
##
## [[952]]
## [1] 1
##
## [[953]]
## [1] 1
##
## [[954]]
## [1] 1
##
## [[955]]
## [1] 1
##
## [[956]]
## [1] 1
##
## [[957]]
## [1] 1
##
## [[958]]
## [1] 1
##
## [[959]]
## [1] 1
##
## [[960]]
## [1] 1
##
## [[961]]
## [1] 1
##
## [[962]]
## [1] 1
##
## [[963]]
## [1] 1
##
## [[964]]
## [1] 1
##
## [[965]]
## [1] 1
##
## [[966]]
## [1] 1
##
## [[967]]
## [1] 1
##
## [[968]]
## [1] 1
##
## [[969]]
## [1] 1
##
## [[970]]
## [1] 1
##
## [[971]]
## [1] 1
##
## [[972]]
## [1] 1
##
## [[973]]
## [1] 1
##
## [[974]]
## [1] 1
##
## [[975]]
## [1] 1
##
## [[976]]
## [1] 1
##
## [[977]]
## [1] 1
##
## [[978]]
## [1] 1
##
## [[979]]
## [1] 1
##
## [[980]]
## [1] 1
##
## [[981]]
## [1] 1
##
## [[982]]
## [1] 1
##
## [[983]]
## [1] 1
##
## [[984]]
## [1] 1
##
## [[985]]
## [1] 1
##
## [[986]]
## [1] 1
##
## [[987]]
## [1] 1
##
## [[988]]
## [1] 1
##
## [[989]]
## [1] 1
##
## [[990]]
## [1] 1
##
## [[991]]
## [1] 1
##
## [[992]]
## [1] 1
##
## [[993]]
## [1] 1
##
## [[994]]
## [1] 1
##
## [[995]]
## [1] 1
##
## [[996]]
## [1] 1
##
## [[997]]
## [1] 1
##
## [[998]]
## [1] 1
##
## [[999]]
## [1] 1
##
## [[1000]]
## [1] 1
##
## [[1001]]
## [1] 1
##
## [[1002]]
## [1] 1
##
## [[1003]]
## [1] 1
##
## [[1004]]
## [1] 1
##
## [[1005]]
## [1] 1
##
## [[1006]]
## [1] 1
##
## [[1007]]
## [1] 1
##
## [[1008]]
## [1] 1
##
## [[1009]]
## [1] 1
##
## [[1010]]
## [1] 1
##
## [[1011]]
## [1] 1
##
## [[1012]]
## [1] 1
##
## [[1013]]
## [1] 1
##
## [[1014]]
## [1] 1
##
## [[1015]]
## [1] 1
##
## [[1016]]
## [1] 1
##
## [[1017]]
## [1] 1
##
## [[1018]]
## [1] 1
##
## [[1019]]
## [1] 1
##
## [[1020]]
## [1] 1
##
## [[1021]]
## [1] 1
##
## [[1022]]
## [1] 1
##
## [[1023]]
## [1] 1
##
## [[1024]]
## [1] 1
##
## [[1025]]
## [1] 1
##
## [[1026]]
## [1] 1
##
## [[1027]]
## [1] 1
##
## [[1028]]
## [1] 1
##
## [[1029]]
## [1] 1
##
## [[1030]]
## [1] 1
##
## [[1031]]
## [1] 1
##
## [[1032]]
## [1] 1
##
## [[1033]]
## [1] 1
##
## [[1034]]
## [1] 1
##
## [[1035]]
## [1] 1
##
## [[1036]]
## [1] 1
##
## [[1037]]
## [1] 1
##
## [[1038]]
## [1] 1
##
## [[1039]]
## [1] 1
##
## [[1040]]
## [1] 1
##
## [[1041]]
## [1] 1
##
## [[1042]]
## [1] 1
##
## [[1043]]
## [1] 1
##
## [[1044]]
## [1] 1
##
## [[1045]]
## [1] 1
##
## [[1046]]
## [1] 1
##
## [[1047]]
## [1] 1
##
## [[1048]]
## [1] 1
##
## [[1049]]
## [1] 1
##
## [[1050]]
## [1] 1
##
## [[1051]]
## [1] 1
##
## [[1052]]
## [1] 1
##
## [[1053]]
## [1] 1
##
## [[1054]]
## [1] 1
##
## [[1055]]
## [1] 1
##
## [[1056]]
## [1] 1
##
## [[1057]]
## [1] 1
##
## [[1058]]
## [1] 1
##
## [[1059]]
## [1] 1
##
## [[1060]]
## [1] 1
##
## [[1061]]
## [1] 1
##
## [[1062]]
## [1] 1
##
## [[1063]]
## [1] 1
##
## [[1064]]
## [1] 1
##
## [[1065]]
## [1] 1
##
## [[1066]]
## [1] 1
##
## [[1067]]
## [1] 1
##
## [[1068]]
## [1] 1
##
## [[1069]]
## [1] 1
##
## [[1070]]
## [1] 1
##
## [[1071]]
## [1] 1
##
## [[1072]]
## [1] 1
##
## [[1073]]
## [1] 1
##
## [[1074]]
## [1] 1
##
## [[1075]]
## [1] 1
##
## [[1076]]
## [1] 1
##
## [[1077]]
## [1] 1
##
## [[1078]]
## [1] 1
##
## [[1079]]
## [1] 1
##
## [[1080]]
## [1] 1
##
## [[1081]]
## [1] 1
##
## [[1082]]
## [1] 1
##
## [[1083]]
## [1] 1
##
## [[1084]]
## [1] 1
##
## [[1085]]
## [1] 1
##
## [[1086]]
## [1] 1
##
## [[1087]]
## [1] 1
##
## [[1088]]
## [1] 1
##
## [[1089]]
## [1] 1
##
## [[1090]]
## [1] 1
##
## [[1091]]
## [1] 1
##
## [[1092]]
## [1] 1
##
## [[1093]]
## [1] 1
##
## [[1094]]
## [1] 1
##
## [[1095]]
## [1] 1
##
## [[1096]]
## [1] 1
##
## [[1097]]
## [1] 1
##
## [[1098]]
## [1] 1
##
## [[1099]]
## [1] 1
##
## [[1100]]
## [1] 1
##
## [[1101]]
## [1] 1
##
## [[1102]]
## [1] 1
##
## [[1103]]
## [1] 1
##
## [[1104]]
## [1] 1
##
## [[1105]]
## [1] 1
##
## [[1106]]
## [1] 1
##
## [[1107]]
## [1] 1
##
## [[1108]]
## [1] 1
##
## [[1109]]
## [1] 1
##
## [[1110]]
## [1] 1
##
## [[1111]]
## [1] 1
##
## [[1112]]
## [1] 1
##
## [[1113]]
## [1] 1
##
## [[1114]]
## [1] 1
##
## [[1115]]
## [1] 1
##
## [[1116]]
## [1] 1
##
## [[1117]]
## [1] 1
##
## [[1118]]
## [1] 1
##
## [[1119]]
## [1] 1
##
## [[1120]]
## [1] 1
##
## [[1121]]
## [1] 1
##
## [[1122]]
## [1] 1
##
## [[1123]]
## [1] 1
##
## [[1124]]
## [1] 1
##
## [[1125]]
## [1] 1
##
## [[1126]]
## [1] 1
##
## [[1127]]
## [1] 1
##
## [[1128]]
## [1] 1
##
## [[1129]]
## [1] 1
##
## [[1130]]
## [1] 1
##
## [[1131]]
## [1] 1
##
## [[1132]]
## [1] 1
##
## [[1133]]
## [1] 1
##
## [[1134]]
## [1] 1
##
## [[1135]]
## [1] 1
##
## [[1136]]
## [1] 1
##
## [[1137]]
## [1] 1
##
## [[1138]]
## [1] 1
##
## [[1139]]
## [1] 1
##
## [[1140]]
## [1] 1
##
## [[1141]]
## [1] 1
##
## [[1142]]
## [1] 1
##
## [[1143]]
## [1] 1
##
## [[1144]]
## [1] 1
##
## [[1145]]
## [1] 1
##
## [[1146]]
## [1] 1
##
## [[1147]]
## [1] 1
##
## [[1148]]
## [1] 1
##
## [[1149]]
## [1] 1
##
## [[1150]]
## [1] 1
##
## [[1151]]
## [1] 1
##
## [[1152]]
## [1] 1
##
## [[1153]]
## [1] 1
##
## [[1154]]
## [1] 1
##
## [[1155]]
## [1] 1
##
## [[1156]]
## [1] 1
##
## [[1157]]
## [1] 1
##
## [[1158]]
## [1] 1
##
## [[1159]]
## [1] 1
##
## [[1160]]
## [1] 1
##
## [[1161]]
## [1] 1
##
## [[1162]]
## [1] 1
##
## [[1163]]
## [1] 1
##
## [[1164]]
## [1] 1
##
## [[1165]]
## [1] 1
##
## [[1166]]
## [1] 1
##
## [[1167]]
## [1] 1
##
## [[1168]]
## [1] 1
##
## [[1169]]
## [1] 1
##
## [[1170]]
## [1] 1
##
## [[1171]]
## [1] 1
##
## [[1172]]
## [1] 1
##
## [[1173]]
## [1] 1
##
## [[1174]]
## [1] 1
##
## [[1175]]
## [1] 1
##
## [[1176]]
## [1] 1
##
## [[1177]]
## [1] 1
##
## [[1178]]
## [1] 1
##
## [[1179]]
## [1] 1
##
## [[1180]]
## [1] 1
##
## [[1181]]
## [1] 1
##
## [[1182]]
## [1] 1
##
## [[1183]]
## [1] 1
##
## [[1184]]
## [1] 1
##
## [[1185]]
## [1] 1
##
## [[1186]]
## [1] 1
##
## [[1187]]
## [1] 1
##
## [[1188]]
## [1] 1
##
## [[1189]]
## [1] 1
##
## [[1190]]
## [1] 1
##
## [[1191]]
## [1] 1
##
## [[1192]]
## [1] 1
##
## [[1193]]
## [1] 1
##
## [[1194]]
## [1] 1
##
## [[1195]]
## [1] 1
##
## [[1196]]
## [1] 1
##
## [[1197]]
## [1] 1
##
## [[1198]]
## [1] 1
##
## [[1199]]
## [1] 1
##
## [[1200]]
## [1] 1
##
## [[1201]]
## [1] 1
##
## [[1202]]
## [1] 1
##
## [[1203]]
## [1] 1
##
## [[1204]]
## [1] 1
##
## [[1205]]
## [1] 1
##
## [[1206]]
## [1] 1
##
## [[1207]]
## [1] 1
##
## [[1208]]
## [1] 1
##
## [[1209]]
## [1] 1
##
## [[1210]]
## [1] 1
##
## [[1211]]
## [1] 1
##
## [[1212]]
## [1] 1
##
## [[1213]]
## [1] 1
##
## [[1214]]
## [1] 1
##
## [[1215]]
## [1] 1
##
## [[1216]]
## [1] 1
##
## [[1217]]
## [1] 1
##
## [[1218]]
## [1] 1
##
## [[1219]]
## [1] 1
##
## [[1220]]
## [1] 1
##
## [[1221]]
## [1] 1
##
## [[1222]]
## [1] 1
##
## [[1223]]
## [1] 1
##
## [[1224]]
## [1] 1
##
## [[1225]]
## [1] 1
##
## [[1226]]
## [1] 1
##
## [[1227]]
## [1] 1
##
## [[1228]]
## [1] 1
##
## [[1229]]
## [1] 1
##
## [[1230]]
## [1] 1
##
## [[1231]]
## [1] 1
##
## [[1232]]
## [1] 1
##
## [[1233]]
## [1] 1
##
## [[1234]]
## [1] 1
##
## [[1235]]
## [1] 1
##
## [[1236]]
## [1] 1
##
## [[1237]]
## [1] 1
##
## [[1238]]
## [1] 1
##
## [[1239]]
## [1] 1
##
## [[1240]]
## [1] 1
##
## [[1241]]
## [1] 1
##
## [[1242]]
## [1] 1
##
## [[1243]]
## [1] 1
##
## [[1244]]
## [1] 1
##
## [[1245]]
## [1] 1
##
## [[1246]]
## [1] 1
##
## [[1247]]
## [1] 1
##
## [[1248]]
## [1] 1
##
## [[1249]]
## [1] 1
##
## [[1250]]
## [1] 1
##
## [[1251]]
## [1] 1
##
## [[1252]]
## [1] 1
##
## [[1253]]
## [1] 1
##
## [[1254]]
## [1] 1
##
## [[1255]]
## [1] 1
##
## [[1256]]
## [1] 1
##
## [[1257]]
## [1] 1
##
## [[1258]]
## [1] 1
##
## [[1259]]
## [1] 1
##
## [[1260]]
## [1] 1
##
## [[1261]]
## [1] 1
##
## [[1262]]
## [1] 1
##
## [[1263]]
## [1] 1
##
## [[1264]]
## [1] 1
##
## [[1265]]
## [1] 1
##
## [[1266]]
## [1] 1
##
## [[1267]]
## [1] 1
##
## [[1268]]
## [1] 1
##
## [[1269]]
## [1] 1
##
## [[1270]]
## [1] 1
##
## [[1271]]
## [1] 1
##
## [[1272]]
## [1] 1
##
## [[1273]]
## [1] 1
##
## [[1274]]
## [1] 1
##
## [[1275]]
## [1] 1
##
## [[1276]]
## [1] 1
##
## [[1277]]
## [1] 1
##
## [[1278]]
## [1] 1
##
## [[1279]]
## [1] 1
##
## [[1280]]
## [1] 1
##
## [[1281]]
## [1] 1
##
## [[1282]]
## [1] 1
##
## [[1283]]
## [1] 1
##
## [[1284]]
## [1] 1
##
## [[1285]]
## [1] 1
##
## [[1286]]
## [1] 1
##
## [[1287]]
## [1] 1
##
## [[1288]]
## [1] 1
##
## [[1289]]
## [1] 1
##
## [[1290]]
## [1] 1
##
## [[1291]]
## [1] 1
##
## [[1292]]
## [1] 1
##
## [[1293]]
## [1] 1
##
## [[1294]]
## [1] 1
##
## [[1295]]
## [1] 1
##
## [[1296]]
## [1] 1
##
## [[1297]]
## [1] 1
##
## [[1298]]
## [1] 1
##
## [[1299]]
## [1] 1
##
## [[1300]]
## [1] 1
##
## [[1301]]
## [1] 1
##
## [[1302]]
## [1] 1
##
## [[1303]]
## [1] 1
##
## [[1304]]
## [1] 1
##
## [[1305]]
## [1] 1
##
## [[1306]]
## [1] 1
##
## [[1307]]
## [1] 1
##
## [[1308]]
## [1] 1
##
## [[1309]]
## [1] 1
##
## [[1310]]
## [1] 1
##
## [[1311]]
## [1] 1
##
## [[1312]]
## [1] 1
##
## [[1313]]
## [1] 1
##
## [[1314]]
## [1] 1
##
## [[1315]]
## [1] 1
##
## [[1316]]
## [1] 1
##
## [[1317]]
## [1] 1
##
## [[1318]]
## [1] 1
##
## [[1319]]
## [1] 1
##
## [[1320]]
## [1] 1
##
## [[1321]]
## [1] 1
##
## [[1322]]
## [1] 1
##
## [[1323]]
## [1] 1
##
## [[1324]]
## [1] 1
##
## [[1325]]
## [1] 1
##
## [[1326]]
## [1] 1
##
## [[1327]]
## [1] 1
##
## [[1328]]
## [1] 1
##
## [[1329]]
## [1] 1
##
## [[1330]]
## [1] 1
##
## [[1331]]
## [1] 1
##
## [[1332]]
## [1] 1
##
## [[1333]]
## [1] 1
##
## [[1334]]
## [1] 1
##
## [[1335]]
## [1] 1
##
## [[1336]]
## [1] 1
##
## [[1337]]
## [1] 1
##
## [[1338]]
## [1] 1
##
## [[1339]]
## [1] 1
##
## [[1340]]
## [1] 1
##
## [[1341]]
## [1] 1
##
## [[1342]]
## [1] 1
##
## [[1343]]
## [1] 1
##
## [[1344]]
## [1] 1
##
## [[1345]]
## [1] 1
##
## [[1346]]
## [1] 1
##
## [[1347]]
## [1] 1
##
## [[1348]]
## [1] 1
##
## [[1349]]
## [1] 1
##
## [[1350]]
## [1] 1
##
## [[1351]]
## [1] 1
##
## [[1352]]
## [1] 1
##
## [[1353]]
## [1] 1
##
## [[1354]]
## [1] 1
##
## [[1355]]
## [1] 1
##
## [[1356]]
## [1] 1
##
## [[1357]]
## [1] 1
##
## [[1358]]
## [1] 1
##
## [[1359]]
## [1] 1
##
## [[1360]]
## [1] 1
##
## [[1361]]
## [1] 1
##
## [[1362]]
## [1] 1
##
## [[1363]]
## [1] 1
##
## [[1364]]
## [1] 1
##
## [[1365]]
## [1] 1
##
## [[1366]]
## [1] 1
##
## [[1367]]
## [1] 1
##
## [[1368]]
## [1] 1
##
## [[1369]]
## [1] 1
##
## [[1370]]
## [1] 1
##
## [[1371]]
## [1] 1
##
## [[1372]]
## [1] 1
##
## [[1373]]
## [1] 1
##
## [[1374]]
## [1] 1
temp2freq(rawd1$MMXT)
## [1] 4760 3990 5770 6270 7660 8910 8960 8600 7630 6520 4290
## [12] 3960 5310 6050 7610 7310 9450 9130 6760 7030 5840 4300
## [23] 3820 4720 4860 6060 7310 7820 8490 8680 7630 6770 5040
## [34] 4350 4590 3810 4840 5850 6910 8150 8210 8720 7340 7080
## [45] 4670 4410 4180 4350 4550 6370 6920 7870 8650 8870 7990
## [56] 6900 6130 4490 3810 3730 5540 6160 6510 7450 8910 8130
## [67] 5930 5210 4180 4540 4570 4060 6430 6820 7460 9020 7930
## [78] 6710 4730 4700 2990 4890 5150 5980 6610 7300 8380 8210
## [89] 7500 7110 5470 4610 4380 4450 4950 6560 6530 7390 9070
## [100] 8430 7940 5980 5420 4040 3710 4490 5040 5330 6650 8010
## [111] 8730 8910 7570 6620 5370 3170 3750 3520 6330 7080 7220
## [122] 8390 9150 8540 7200 6890 5140 4390 4180 3770 5870 5950
## [133] 6900 8130 8600 8230 7190 5960 4550 3990 4220 4610 3820
## [144] 6030 6920 8150 8210 8130 6770 5800 5540 4290 4090 3360
## [155] 4430 6370 6990 7950 8260 8900 7760 5750 5480 4070 4750
## [166] 4200 5470 6060 7400 7730 9010 8750 7690 6420 6030 3610
## [177] 4110 4810 5270 6980 6640 7280 7970 8920 6930 7000 4970
## [188] 4160 2190 4380 5630 6070 6360 7590 8980 8660 7730 5930
## [199] 4600 3580 3630 3870 4200 5410 6770 7850 9490 8880 8060
## [210] 6530 6060 4880 3790 4010 5570 5640 6980 8730 8580 8450
## [221] 7550 6860 4800 4540 4880 3900 5150 6650 7540 9030 9720
## [232] 9170 7710 5260 4400 3420 4440 4270 4500 5220 6820 7900
## [243] 8960 8820 7880 6120 4640 4440 4530 4930 5450 5950 6820
## [254] 8410 9070 8940 7010 7150 5120 3840 3330 3170 4820 5760
## [265] 6950 8500 8680 9160 8390 7080 4580 3810 4470 3880 5080
## [276] 6550 7220 7840 9270 8870 7960 5900 5700 4140 3790 4880
## [287] 4350 5710 7510 7940 8810 8880 7840 6450 4890 3250 4390
## [298] 5140 5340 6410 7480 7920 9160 8810 7280 5320 4830 4640
## [309] 4220 4750 5150 6760 7190 8170 9090 8760 6780 6800 5180
## [320] 3910 3860 4360 5180 5830 6260 7960 8890 8260 7450 6970
## [331] 5130 3290 3990 4110 5470 6150 7960 7530 8960 8790 8030
## [342] 6390 5240 4130 3030 3530 5070 5910 7120 8070 9260 9470
## [353] 7280 6780 4870 4840 2110 5090 4880 7240 7720 8460 9140
## [364] 8920 8010 5880 5290 4120 4480 4570 5080 6510 7880 8540
## [375] 9340 9250 8080 6820 4680 4180 3560 4400 4260 6410 7510
## [386] 8140 9140 8960 8140 6230 5790 3430 4260 3580 5420 6050
## [397] 6940 8800 9590 8650 7900 7200 5700 5280 5190 5100 6060
## [408] 7300 8360 8150 9450 9360 7610 7090 5830 4430 4220 4820
## [419] 4810 5840 6790 8060 9270 9130 8360 6890 4340 4280 3980
## [430] 2520 5200 6450 8150 8420 9690 9130 8090 6880 4510 4610
## [441] 2030 4140 4490 5920 7660 7950 9120 9150 8160 7160 5300
## [452] 4460 4370 4570 5140 6190 6830 8180 8750 8910 8790 6510
## [463] 4500 4680 4620 3700 5340 6910 7770 7580 9100 9060 8060
## [474] 6740 6250 5060 3380 4370 5520 6190 7830 8450 9230 9240
## [485] 8010 7050 4500 4740 4390 4680 5660 6300 7600 7990 8960
## [496] 8730 7090 6480 5530 4220 3520 3540 4750 6440 6760 7580
## [507] 9010 8920 7900 6840 5030 4600 3750 4680 4510 6870 6710
## [518] 7610 8930 8930 8210 6970 5470 4350 4200 4330 4710 6490
## [529] 7600 7570 8640 8730 7970 7390 5150 4200 4530 4670 5410
## [540] 5540 6990 7370 9070 9040 7550 7260 4920 3920 4210 4630
## [551] 5630 7190 6830 8150 9350 9000 7900 5960 5060 4940 4230
## [562] 4810 5410 6500 7790 7650 9540 9210 7860 7130 4580 4530
## [573] 4380 4160 4710 6460 7420 8090 9060 7150 4890 3780 2760
## [584] 4210 5010 7200 7630 7920 9210 9430 8020 6040 6430 4550
## [595] 3570 5080 4960 6240 3480 4370 4560 6410 7420 7620 9220
## [606] 8630 7580 6170 4970 3940 4010 4420 4620 6950 7250 8320
## [617] 8980 9100 8450 7370 4840 4510 5250 4650 5670 5800 6750
## [628] 8040 9400 9300 8520 7270 5840 4720 4100 5400 4620 6620
## [639] 7410 7710 9450 8850 8040 6690 6160 4550 3790 3960 4380
## [650] 5740 7250 7810 8820 9440 7820 7050 4270 4190 4140 4050
## [661] 5120 6150 7240 8630 9090 8620 8380 6890 5150 4940 3080
## [672] 4620 5250 6120 7380 8110 9180 9160 8110 6400 4840 5140
## [683] 4680 5030 4770 6050 8380 8070 8630 9250 8070 7270 5350
## [694] 4680 4540 4160 5310 6470 6760 8450 9250 9060 7640 6510
## [705] 5030 4820 3990 4040 5330 6250 7380 8420 9650 9010 8330
## [716] 7000 5280 4560 4840 5550 5840 6130 7470 9110 9600 9610
## [727] 7050 6500 5030 4160 3670 4380 5010 6980 7140 8300 8760
## [738] 8870 7900 7080 5580 5000 3280 5410 5690 6300 7420 8040
## [749] 9290 9290 8530 7520 5580 4280 4390 4690 4760 6210 7280
## [760] 7910 9600 8860 7850 7360 5200 4520 4830 4890 4440 6430
## [771] 7020 8100 9100 8650 6540 7580 5690 4830 4510 4560 5600
## [782] 6260 7980 7980 9430 9010 8480 6810 5490 4570 4830 5010
## [793] 5000 5830 7240 7780 9280 9420 8760 6910 5350 3870 4380
## [804] 4830 6020 6070 6980 8040 9200 8450 7910 6970 5180 4040
## [815] 3530 4320 4830 6990 7810 7860 9050 9410 8320 5800 5710
## [826] 4810 4310 5430 4920 5520 7440 8580 9290 9550 7530 6420
## [837] 5340 4350 4240 4910 5290 6270 7420 8160 9030 9730 7680
## [848] 6560 5340 4060 3910 5020 6290 6360 7290 8560 8860 9300
## [859] 7640 6500 5120 3670 3950 4590 5410 5920 7460 8340 9310
## [870] 9280 7650 7050 4940 4790 3980 5110 5320 6750 7040 8960
## [881] 9690 8710 8170 7290 5470 4450 4190 4190 4860 5200 6860
## [892] 7890 9320 8630 8110 6650 5140 4950 4630 5010 5070 6470
## [903] 7930 8240 9280 8890 8240 6910 5600 4970 4030 5380 5040
## [914] 7160 7250 8630 9030 8830 7960 7000 5080 4510 4160 4550
## [925] 5840 6510 7060 8350 8940 8990 8200 7230 4410 3680 2720
## [936] 4530 5310 6290 7400 8430 9330 9070 8990 7300 4840 5310
## [947] 3770 4800 4990 7010 7610 8340 9150 8650 8020 6980 5480
## [958] 5430 5170 4860 5990 6680 7120 8000 9280 9460 8440 6520
## [969] 5770 4610 4000 4360 5300 5990 7070 8000 8900 9370 7760
## [980] 6560 4750 4250 5120 5340 5810 6170 7190 8190 9040 9440
## [991] 7820 6980 5260 2890 4550 5130 5430 6330 7500 8060 9200
## [1002] 9400 7580 6250 5370 4120 3800 4290 5170 7020 7920 8640
## [1013] 9730 8890 7200 6620 3760 4050 5100 4690 6500 6550 7540
## [1024] 8940 8740 9260 7150 7030 5080 4510 4530 5180 5920 7690
## [1035] 7940 8690 8790 8690 8540 7570 5580 4490 4190 5130 5610
## [1046] 6990 7910 9360 9750 9480 8080 7620 5110 4390 4620 3110
## [1057] 5240 6560 7390 8500 9620 8870 8240 6720 5660 4430 4730
## [1068] 4940 5560 6880 7220 8480 9340 9210 8870 7020 5840 3680
## [1079] 4290 5800 5530 6210 7120 8200 9490 9700 8350 6830 4850
## [1090] 4950 4970 5700 6560 7060 8070 8500 8540 8950 8210 7320
## [1101] 4970 3840 3900 4010 6000 6510 7910 7970 8020 8450 8000
## [1112] 6940 4850 4810 5140 4690 6250 6970 7950 8550 9330 9550
## [1123] 8730 6580 4830 4650 4450 5280 5500 6310 7150 8320 9210
## [1134] 9550 5680 4770 4040 4760 5050 6740 6770 8680 9440 9460
## [1145] 8000 6780 4950 4350 4080 4660 5870 5880 7690 8170 8780
## [1156] 8960 8350 6730 5300 4280 4490 5060 5430 6650 7910 7620
## [1167] 9410 9410 8720 6780 5580 4260 4720 5210 5910 6300 7200
## [1178] 8040 9230 9420 7900 7420 6510 4870 4660 5410 5990 7050
## [1189] 7480 8420 9600 9540 8060 6700 4220 4210 4430 4260 5710
## [1200] 6550 8060 8450 9310 9780 8790 6870 6190 4580 4590 4820
## [1211] 4900 6550 7300 8380 9650 8960 8280 6440 5510 5050 5140
## [1222] 4450 5640 6810 7480 8360 9990 9860 8430 7690 4950 4860
## [1233] 4470 4830 6490 7080 7220 8120 9240 8960 7930 6940 5560
## [1244] 5050 4540 5580 5970 6790 7410 8140 9650 9240 8360 7340
## [1255] 5420 4260 5210 4680 5620 7110 7880 8740 9910 9250 8160
## [1266] 6770 5450 4760 4350 5030 6450 6730 7940 8780 10230 9450
## [1277] 8270 6950 5530 4350 4360 5070 5450 6380 7220 8270 9500
## [1288] 9440 8150 7110 5860 3980 4790 5250 5370 6360 7830 8090
## [1299] 9250 9120 9050 6030 5650 3700 4400 4490 5820 6400 6690
## [1310] 7820 8980 8850 8110 7280 4700 4250 4340 4080 5440 5950
## [1321] 6840 7940 9340 9300 8650 6870 5230 4620 5050 4620 6370
## [1332] 7030 7260 8510 9610 9430 8560 6700 5690 4400 4390 4710
## [1343] 5590 6180 7380 8420 9420 9300 8060 6300 5390 3990 4840
## [1354] 3690 5230 6500 7400 7830 9420 8900 8250 7620 5050 4800
## [1365] 5490 6480 6710 7090 8980 9010 9200 8390 7280 4860
for (i in 1:164) {
play_short(temp2freq(rawd1$MMXT[i]))
}